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I. INTRODUCTION 

Electrical current is the transfer of charge from one 
end of the conductor to the other. The statistics of this 
charge transfer was investigated by Levitov and Lesovik 
jjj. It is binomial for a single-channel conductor at zero 
temperature and double-Poissonian at finite temperature 
in the tunneling regime [||. The second cumulant, the 
noise power, has been measured in a variety of systems 
H . Ways of measuring the third cumulant have been 
proposed |4|], but not yet carried out. 

Electrical current also transfers momentum to the lat- 
tice. The second cumulant, the electromechanical noise 
power, determines the mean square displacement of an 
oscillator through which a current is driven. It has been 
studied theoretically @, ^, ^, || and is expected to lie 
within the range of sensitivity of nanomechanical oscil- 
lators ||. No theory exists for higher order cumulants 
of the transferred momentum (which would determine 
higher cumulants of the oscillator displacement). It is 
the purpose of the present paper to provide such a the- 
ory. 

In the context of charge transfer statistics there exist 
two approaches: a fully quantum mechanical approach 
using Keldysh Green functions jl], |l0j and a semiclassical 
approach using the Boltzmann-Langevin equation 
Here we take the former approach, to arrive at a quantum 
theory of momentum transfer statistics. As a test, we 
show that the second moment calculated from Keldysh 
Green functions coincides in the semiclassical limit with 
the result obtained from the Boltzmann-Langevin equa- 
tion by Shytov, Levitov, and one of the authors ||. 

A calculation of the complete cumulant generating 
function of transferred momentum (or, equivalently, of 
oscillator displacement) is presented for the case of a 
single-channel conductor with a localized scatterer. The 
generating function in this case can be written entirely 
in terms of the transmission probability P of the scat- 
terer. In the more general multi-channel case one also 
needs knowledge of the wavefunctions. This is an essen- 
tial difference from the charge transfer problem, which 
can be solved in terms of transmission eigenvalues for 
any number of channels. At zero temperature the mo- 
mentum statistics is binomial, just as for the charge. At 



finite temperature it is multinomial, even in the limit 
r — > 0, different from the double-Poissonian distribution 
of charge. 

The outline of the paper is as follows. In Sec. II we 
formulate the problem in a way that is suitable for fur- 
ther analysis. The key technical step in that section is a 
unitary transformation which eliminates the dependence 
of the electron-phonon coupling Hamiltonian on the (un- 
known) scattering potential of the disordered lattice. The 
resulting coupling Hamiltonian contains the electron mo- 
mentum flow and the phonon displacement. In the next 
section we use that Hamiltonian to derive a general for- 
mula for the generating function of the distribution of 
momentum transferred to a phonon (as well as the dis- 
tribution of phonon displacements) . It is the analogue of 
the Levitov-Lesovik formula for the charge transfer distri- 
bution 0| . For a localized scatterer we can evaluate this 
statistics in terms of the scattering matrix. We show how 
to do this in Sec. IV, and give an application to a single- 
channel conductor in Sec. [y| In Sees. VI and VII we turn 
to the case that the scattering region extends throughout 
the conductor. We follow the Keldysh approach to derive 
a general formula for the generating function and check 
its validity b y red eriving the result of Ref . M . We con- 
clude in Sec. VIII with an order of magnitude estimate of 



higher order cumulants of the momentum transfer statis- 
tics. 



II. FORMULATION OF THE PROBLEM 

The excitation of a phonon mode by conduction elec- 
trons is described by the Hamiltonian 

h = fiata + pV 2m + v ^ - <2 u ( r *)b 

i i 

where we have set fi = 1. The phonon mode has anni- 
hilation operator a, frequency fl, mass M, and displace- 
ment Qu(r), where Q = (2Mf2) -1 / 2 (a + a^) is the ampli- 
tude operator. The electrons have position rj, momen- 
tum p s ; = —id/dri, and mass m. Electrons and phonons 
are coupled through the ion potential V(r). We assume 
zero magnetic field. Electron-electron interactions and 
the interactions of electrons and phonons with an exter- 
nal electric field have also been omitted. 
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We assume that electrons and phonons are uncoupled 
at time zero and measure moments of the observable 
A of the phonons after they have been coupled to the 
electrons for a time t. The operator A(a, cv) could be 
the amplitude Q of the phonon mode, its momentum 
a — aJ), or its energy fla^a. The mo- 
ment generating function for A is 



m=0 



tin 

m = Y,h {Am{t)) = TieiAe ~ 



lHt pe iHt . 



(2.2) 



The initial density matrix p — p e p P is assumed to factor- 
ize into an electron and a phonon part. 

We assume small displacements, so an obvious way to 
proceed would be to linearize V(r — Qu) with respect to 
the phonon amplitude Q. Such a procedure is compli- 
cated by the fact that the resulting coupling — Qu • W 
of electrons and phonons depends on the ion potential V . 
Because of momentum conservation, it should be possi- 
ble to find the momentum transferred by the electrons 
to the lattice without having to consider explicitly the 
force — W. In the semiclassical calculation of Ref. ||] 
that goal is achieved by the continuity equation for the 
flow of electron momentum. The unitary transformation 
that we now discuss achieves the same purpose in a fully 
quantum mechanical framework. 

What we need is a unitary operator U such that 



U^V[r - Qu(r)]U = V(r). 



(2.3) 



For constant u we have simply U = cxp[— iQu ■ p]. More 
generally, for space-dependent u, we need to specify the 
operator ordering (denoted by colons : • • • :) that all po- 
sition operators r stand to the left of the momentum 
operators p. We also need to include a Jacobian deter- 
minant || J|| to ensure unitarity of U. As shown in App. 
[A|, the desired operator is 



U = 



1/2 . e -iQu(r)-p 



J a p = <W - Qd a up(r), 



(2.4) 



with d a = d/dr a . All this was for a single electronic 
degree of freedom. The corresponding operator for many 
electrons is U = Ili^> where Ui is given by Eq. (2.4) 
with r, p replaced by r i? p,. 

The Hamiltonian ( |2.l| ) transforms as HU 
H m t, with 



H 



H = Qa^a + ^[p?/2m + V(r 2 )], 

i 

i 



H 



-QF-^jPTl + 0(u 2 ). 



(2.5a) 
(2.5b) 



Here F is the driving force of the phonon mode, 

F = J— y][u a /3(Ti)piaPii3 + Pi a U a p(ri)pi0] + H.C, 



and II is the total electron momentum, 

n=|5^u(r i ).p i +H.c, (2.7) 



weighted with the (dimensionless) mode profile u(r). We 
have defined the shear tensor u a p — h(d a up + dpu a ). 
The abbreviation H.c. indicates the Hermitian conjugate 
and a summation over repeated cartesian indices a, (3 is 
implied. 

The interaction Hamiltonian H m t is now independent 
of the ion potential, as desired. In the first term — QF 
we recognize the momentum flux tensor, while the second 
term — PH/M is an inertial contribution to the momen- 
tum transfer. The inertial contribution is of relative or- 
der flX/vF (A being the wavelength of the phonon and vf 
the Fermi velocity of the electrons) and typically <C 1 . In 
what follows we will neglect it. We also neglect the terms 
in Hi n t of second and higher order in u, which contribute 
to order Xf/L to the generating function (with L the 
length scale on which u varies). These higher order in- 
teraction terms account for the momentum uncertainty of 
an electron upon a position measurement by the phonon. 

If we apply th e un itary transformation U to the gen- 
crating function (2.2), we need to transform not only H 
but also A — > U^AU = A and p — ► U'pU = p, resulting 



^(0 - Tre^ A e-^ Ho+H ^pe^ Ho+H - 



(2.8) 



In App. |X| we show that, quite generally, the distinction 
between /?, A and p, A is irrelevant in the limit of a long 
detection time t, and we will therefore ignore this dis- 
tinction in what follows. 

If u is smooth on the scale of Xf, so that gradients of 
u a /3 can be neglected, one can apply the effective mass 
approximation to the Hamiltonian (2.5). The ion poten- 
tial V = Vi a t + Vi mp is decomposed into a contribution 
Vi a t from the periodic lattice and a contribution Vi mp 
from impurities and boundaries that break the periodic- 
ity. The effects of Vj a t can be incorporated in an effec- 
tive mass m* (assumed to be deformation independent 
|fl2| , p^[ ) and a corresponding quasimomentum p*. The 
unperturbed Hamiltonian takes the usual form 

H = r! a t a + ^[ p *2/2m* + ^ mp (r l )]. (2.9) 



As shown in App. [B], the force operator in Hint is then 
expressed through the flow of quasi-momcntum, 



(2.10) 



(2.6) 



whereas the inertial contribution is still given by Eq. (2.7) 
in terms of the true electron momentum. 
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III. MOMENTUM TRANSFER STATISTICS 



A. Generating function 



The role of the integrated current J(t) = f*dt'I(t r ) is 
taken in our problem by the integrated force K(t). 



A massive phonon mode absorbs the momentum that 
electrons transfer to it without changing its displacement. 
We may therefore define a statistics of momentum trans- 
fer to the phonons without back action on the electrons 
by choosing the observable A = P = — i(MQ/2) 1 ' 2 (a — 
a^) in Eq. ( |2.2| ) and taking the limit M — > oo, £1 — > at 
fixed MCI. We assume that the phonon mode is initially 
in the ground state, so that ap p = 0. 

We transform to the interaction picture by means of 
the identity 



iH t-iHt 



= Texp 



dt'H int (t r ) 



(3.1) 



where T denotes time ordering (earlier times to the right 
of later times) of the time dependent operator Hi nt (t) = 
e zHot Hi nt e~ lHot . In the massive phonon limit we have 
-ffint(i) = —QF(t) with time independent Q (since Q 
commutes with Hq when O — > 0) . Eq. (2.2) takes the 
form 



B. Relation to displacement statistics 

Cumulants ((AP(i))) of the momentum transferred 
in a time t are obtained from the cumulant generat- 
ing function ln^(0 = E„(( AP ( t )"»^"/ n! - Cumu- 
lants ((F{ui) n )) of the Fourier transformed force F(u) — 
Jdte luJt F(t) then follow from the relation AP(t) = 

J Q dt' F(t'). The limit t — > oo of a long detection time 
corresponds to the low frequency limit, 

//n^)^ - 2tt* f$>^ Inn- «AP(t)™)) . 

(3.6) 

Cumulants of the Fourier transformed displacement 
Q(w) of the oscillator follow from the phenomenological 
equation of motion 



HO = (T ± ex.p[-iQK_(t)]e SP ex V [iQK + (t)]), (3.2) Q(w) = R(lo)F{u), R(u) 



1 

M 



[ijj 2 — U! 2 — ILO(jJ c /Q)~ 



where K±(t) = f dt± F(t±) and 7± denotes the Keldysh 
time ordering: times t_ to the left of times t + , earlier t_ 
to the left of later t_, earlier t+ to the right of later t+. 

Taking the expectation value of the phonon degree of 
freedom we find 



(3.7) 



_ p a 2 Mn/2 



x ( T± exp 



(K+ - K-) 2 



(3.3) 



Here M. is the motional mass, lo c the characteristic fre- 
quency, and Q the quality factor of the oscillator. Since 
the force noise is white until frequencies that are typically 
3> w c , one has in good approximation 



\i=l / 3 = 1 




The factor exp(£ 2 Mil/2) originates from the uncertainty 
(Affi) 1 / 2 of the momentum of the phonon mode in 
the ground state (vacuum fluctuations). It is a time- 
independent additive contribution to the second cumu- 
lant and we can omit it for long detection times. The 
quadratic term cx K±/MQ, becomes small for a small un- 
certainty (Mil) -1 / 2 of the displacement in the ground 
state. It describes a back action of the phonon mode on 
the electrons that persists in the massive phonon limit. 
(A similar effect is known in the context of charge count- 
ing statistics JT^].) This term may be of importance in 
some situations, but we will not consider it here, assum- 
ing that the electron dynamics is insensitive to the vac- 
uum fluctuations of the phonon mode. 

With these simplifications we arrive at a formula for 
the momentum transfer statistics, 



x lim - «AP(t)"» . (3. 

t— >oo t 



Optical or magnetomotive detection of the vibration, 
as m Refs. |J, llq , 0, measures the probability distri- 
bution P(Q) of the displacement at any given time. The 
cumulants o f P (Q) are obtained by Fourier transforma- 
tion of Eq. (|H), 



{{Q n )) = K n ]im -«AP(t)"», 



(3.9) 



/d f duj f "rt> \ 

27 "J -^ R (^--- R M2n6 (X>0 



dtR(t) 7 



(3.10) 



HO = {T ± ew[&K-(t)]exp[&K+{t)]), 



(3.4) 



that is of the same form as the formula for charge count- 
ing statistics due to Levitov and Lesovik JD], 

charged) = (T ± exp[i£J_(i)]exp[|£J + (t)]). (3.5) 



For Q 3> 1 the odd moments can be neglected, while the 
even moments are given by 



1 



2k uj c 



(3.11) 
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FIG. 1: Sketch of a freely suspended wire. The matrices t, if 
and r, r describe transmission and reflection by a localized 
scatterer (shaded). A voltage V drives a current through the 
conductor, exciting a vibration. 



EVALUATION IN TERMS OF THE 
SCATTERING MATRIX 



IV. 



The Levitov-Lesovik formula ( p.5| ) for the charge trans- 
fer statistics can be evaluated in terms of the scatter- 
ing matrix of the conductor |[], [l8|, without explicit 
knowledge of the scattering states. This is possible be- 
cause the current operator depends only on the asymp- 
totic form of the scatter ing states, far from the scatter- 
ing region. The formula (3.4) for the momentum transfer 
statistics can be evaluated in a similar way, but only if 
the mode profile u(r) is approximately constant over the 
scattering region. 

To this end, we first write the force operator ( |2.6| ) in 
second quantized form using a basis of scattering states 

Vw( r ), 



F(t) = 



dede' 
2tt 



Y / e li£ - £ ' )t ci(e)M nn ,(s,s')c nl (e'), 

(4.1) 



M nn >(e,e') = — I drip* iE (p a u a ppp 

Tfl J 

+ [[Ua0,p a ],Pp])lp n r ie >. (4.2) 



The operator c„(e) annihilates an electron in the n-th 
scattering channel at energy e. The mode index n runs 
from 1 to N (or from N + 1 to 2N) for waves incident 
from the left (or from the right). (Se e Fig. [j] for a dia- 
gram of the geometry and see Ref. [M for the analogous 
representation of the current operator.) The commuta- 
tor [[u aj 3,p a ],P/3] can be neglected if u is smooth on the 
scale of the wavelength (hence if Xp/L -C 1). 

We assume that the derivative u a p of the mode profile 
vanishes in the scattering region, so that for the scatter- 
ing states we may use the asymptotic form 



^ n , e (r) = < e (r) + £ S m „( £ )C u », 



(4.3) 



in terms of incident and outgoing waves (f>lf'° ut (normal- 
ized to unit current) and the scattering matrix S rnn (e). 
Since we are neglecting the Lorentz force we may assume 
that 4)fff\ — <J>m* £ - The scattering matrix has the block 



£OUt 

structure 



S = 



r t' 
t r' 



(4.4) 



with Nx N transmission and reflection matrices t, t', r, r' . 
These matrices are related by unitarity (S — S') and 
possibly also by time-reversal symmetry (S — S T ). 

The operator p a u a ppp will couple only weakly the inci- 
dent to the outgoing waves, provided u is smooth on the 
scale of \f, and we neglect this coupling. The matrix M 
then separates into an incident and outgoing part, 

M(e,e') = M in (e,e') + 5 t (e)M out (e, e')S(e'). (4.5) 



The matrices M m and M out are defined as in Eq. (gj) 
with tp replaced by </> m and 4>° nt , respectively. (They are 
Hermitian and related by M out = M 111 *.) These two ma- 
trices vary with energy on the scale of the Fermi energy 
Ep, while the scattering matrix S has a much stronger 
energy dependence (on the scale of the Thouless energy). 
We may therefore replace M m , M out by their value at 
e = e' = Ep and assume that the energy dependence of 
M is given entirely by the scattering matrix. 

The force operator can similarly be separated into 
F = F in + F out , where F in and F out are defined as in Eq. 
(O) with the matrix M replaced by M in and S f M out S, 



respectively. We now proceed in the same way as in Ref. 
]l9| for the current operator, by noting that the analytic- 
ity of S(e) in the upper half of the complex plane implies 
simple commutation relations: 



[F in (t),F in (t')} = 0, [F out (t),F out 
[F [n (t),F out (t')} =0 if t >t'. 



(*')] = 0, Vt,t', 
(4.6) 



It follows that the Keldysh time ordering T± of the force 
operators is the same as the so-called input-output or- 
dering, defined by moving the operators F ln (t-) to the 
left and F- m (t + ) to the right of all other operators — irre- 
spective of the value of the time arguments. The reason 
for preferring input-output ordering over time ordering 
is that Fourier transformation from time to energy com- 
mutes with the former ordering but not with the latter. 

In the limit t — ► oo different energies become uncou- 
pled, and the cumulant generating function takes the sim- 
ple form 



t 

2^ 



de ln{e Fin ^ 2 e F °^e F " 1 ^ 2 ) 



(4.7) 



entirely analogous to the input-output ordered formula 
for charge transfer . The Fourier transformed force is 
defined as 

F in (e) = c t (e)Af in (£,£)c(e), (4.8a) 
F out (e) = c\s)S^e)M out (e,e)S(e)c(e). (4.8b) 
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(The operators c„ have been collected in a vector c.) 
The matrices M m ' out are block diagonal, 



M m = M c 



M L 
M R 



(4.9) 



but the N x N matrices M^r are in general not di- 
agonal themselves. They take a simple form for a lon- 
gitudinal phonon mode, when u is a function of x in 
the x-direction (along the conductor), so that u Q ^(r) = 
SaxSf3xu'(x). The commutator [[u',pa;],Pa;] does not con- 
tribute because </>™' out is an eigenstate of p x (with eigen- 
value p™ = - 
tion one has 



p° ut = p n ). Hence for a longitudinal vibra- 



(M L )„ 
(M R ) r 



= 5 n 



>\p n \(u - U L ), 
>\p n \(u R ~ Uq). 



(4.10a) 
(4.10b) 



The value of u(x) in the scattering region is denoted by 
Mo, while ul,ur denote the values at the left and right 
end of the conductor. The more complex situation of a 
transverse phonon mode, when the matrices M Lj r are no 
longer diagonal, is treated in Ref . ]2l| . 

W e ar e now ready to calculate the expectation value in 
Eq. (4.7). We assume that the incident waves originate 
from reservoirs in thermal equilibrium at temperature T, 
with a voltage difference V between the left and right 
reservoir. The Fermi function in the left (right) reservoir 
is /l (/r). We collect the Fermi functions in a diagonal 
matrix / and write 

(ct(e)c nl (e'))=f nn ,(e)5(e-s'), f = ^ £ 

(4.11) 

All o the r expectation values of c, c* vanish. We evaluate 
Eq. (4.7) with help of the determinantal identity 



JJexp( c tA lC )) = l-/ + /JJ e 



(4.12) 



valid for an arbitrary set of matrices A4, and the identity 
exp(S* t A5) = SU A S, (4.13) 
valid for unitary S. The result is 

Xnm = h J de\n\\l-f + fei Mia SHe)^ Mmt S(e)\\, 

(4.14) 

where we have also used that the two matrices M ln and 
/ commute. 

At zero temperature /l = 9(EF + eV— e), /r = 9{Ef — 
e). The energy range e < Ep, where /l = /r = 1, 
contributes only to the first moment, while the energy 
range Ep < e < Ep + eV, where / L = 1 and /r = 
0, contributes to all moments. For small voltages we 
may neglect the energy dependence of S(e) in that range. 



Using the block structure @, Q, of S, Af in < out the 
generating function for the second and higher cumulants 
takes the form 

pVt 

ln^) = ^-\ii\\rh^r + th^t\\+O{0. (4.15) 

(By 0(£) we mean terms linear in £.) This determinant 
can not be simplified further without knowledge of S. 
That is a major complication relative to the analogous 
formula for the charge transfer statistics Q, which can 
be cast entirely in terms of the transmission eigenvalues 
T„ (eigenvalues of tt^): 



In -^charge (£) 



2vr 



ds 



+ r n (e e «-l)/ L (l-/ R ) 

+ r n (e- e «-l)/ R (l-/ L )] . (4.16) 

In the case of momentum transfer, eigenvalues and eigen- 
vectors both play a role. 



APPLICATION TO A ONE-DIMENSIONAL 
CONDUCTOR 

A. Straight wire 



Further simplification of Eqs. (4.14) and (4.15) is pos- 



sible if the conductor is so narrow that it supports only 
a single propagating mode to the left and right of the 
scattering region (N = 1). The scattering matrix then 
consists of scalar transmission and reflection coefficients 
t, t', r, r' (related to each other by unitarity). We consider 
the case of a longitudinal vibration with 



M m = M out = pp 



u -ul 

MR - Uq 



(5.1) 



cf. Eq. (4.1C). Because of unitarity the result depends 
only on the transmission probability T = \t\ 2 = \t'\ 2 — 



„'|2 



ln^(0 = ^J d£ H 1 + (e 2 ^" R -" L > - 1)/ L / R 
+ r(e &>H«*-« L ) _ _ /r) + /r(1 _ 

+ (1 - r)(e 2 ^("«-" L ) - 1)/ L (1 - / R ) 

+ (1 - r)(e 2 ^^-"«) - 1)/ R (1 - /l)] • (5.2) 
At zero temperature this simplifies further to 
In = tl l n [l + re «PF(uR+«L-2u ) _ pj + 

(5.3) 

The zero temperature statistics ( |5.3| ) is binomial, just 
as for the charge. [The generating function Charge (£) 
at T — is obtained from Eq. (|5.3D after substitution 



of pf{ur + Ul — 2u ) by e, cf. Eq. ( 4.16 ).] At finite 
temperatures one has the multinomial statistics (5.2), 
made up of stochastically independent elementary pro- 
cesses with more than two possible outcomes. The ele- 
mentary processes may be characterized by the numbers 
(n^,n?;) € {0,1} of electrons incident on the scatterer 



from the left, right and the numbers [tit. 



t)G{0,l} 



of outgoing electrons to the left, right. The non- vanishing 



probabilities P [ (n^ , nf^ ) 
events evaluate to: 



(n, 



L R 
out s "'out 



)] of scattering 



P[(0,0) 

mi) 
mi) 

p[0-,o) 
p[0-,o) 

mi) 



(0,0)] 
(0,1)] 
(1,0)] 

(1.0) ] 
(0,1)] 

(1.1) ] 



(1-/l)(1-/r), 

(i-/L)/ R (i-r), 

(1 - /l)/rT, 

/ L (i-/ R )(i-r), 
/ L (i-/ R )r, 



(5.4) 



Thes e probabilities appear in the generating function 
(5.3), multiplied by exponentials of £ times the amount 
of transferred momentum. 

A longitudinal vibration of a straight wire clamped at 
both ends would corr espo nd to ul = Mr = 0, u g ^ . In 
that special case Eq. (5.2) is equivalent to Eq. ( 4.16 ) for 
-?charge(C) under the substitution r — > 1 — T, 2ppUo — > e. 
In this case the multinomial statistics becomes a double- 
Poissonian in the limit T — ► 0, corresponding to two inde- 
pendent Poisson processes originating from the left and 
right reservoirs 0. A longitudinal vibration is difficult 
to observe, in contrast to a transverse vibration which 
can be observed optically J15L [l6| or magnetomotively 
p7[ . However, the direct excitation of a transverse mode 
is not possible in a single-channel conductor, while in a 
multi-channel conductor (width W) it is smaller than the 
excitation of a longitudinal mode by a factor (W/L) 2 [ fjl|| . 
So it would be desirable to find a way of coupling lon- 
gitudinal electron motion to transverse vibration modes. 
In the following subsection we discuss how this can be 
achieved by bending the wire. 



B. Bent wire 

The bending of the wire is described as explained in 
Ref. p2[ , by means of a vector f2(s) that rotates the lo- 
cal coordinate system e x (s), e y (s), e z (s) as one moves an 
infinitesimal distance ds along the wire: 5e a = Q x e a 5s. 
The local coordinate x is along the wire and y, z are per- 
pendicular to it. The component f2i of Jl along the wire 
describes a torsion (with |f2ii| the torsion angle per unit 
length) , while the perpendicular component £l± describes 
the bending (with ISIjJ 1 the radius of curvature). 

The momentum operators and wavefunctions, written 
in local coordinates, depend on the bending by terms 
of order A^|f2|, which we assume to be <C 1. These 
quantities may therefore be evaluated for a straight wire 
(f2 = 0). The dependence on the bending of the strain 







FIG. 2: Two vibration modes in a bent wire (top) and the 
corresponding longitudinal displacements u c ff in the straight 
wire (bottom). 



tensor is of order L\Ct\ an d ca n not be neglected. For the 
interaction Hamiltonian ( |2.5|) we need Vu in the global 
coordinate system. It is obtained by differentiating the 
local coordinates of u as well as the local basis vectors. 
A bent wire can then be represented by a straight wire 
with an effective displacement u c g related to u (in local 
coordinates) by 



d 


d 


-SI x u, 




TT "off = 

OX 






dx 






d 


d 


d 


d 


oy 


oy 


— Ueff = 

Oz 


Oz 



(5.5a) 
(5.5b) 



The second term on the right-hand-side of Eq. ( |5.5a|) ac- 
counts for the centrifugal force exerted by an electron 
moving along the bent wire. It rotates a transverse mode, 
with u pointing in radial direction, into a fictitious lon- 
gitudinal mode with u c s^ x of order L|f2j_|. 

Fig. H shows two vibration modes in a bent wire with 
the corresponding longitudinal component u c $ iX of the ef- 
fective displacement. To apply the formulas of the previ- 
ous subsection we need ul = u e s >x {xi), ur. = UcS x{xb), 
and uq — u e s tX (xo)- The first mode, Fig. |ja, has 
ul = mr = and % ^ 0. It measures the amount 
of electron momentum that has been transferred to the 
scatterer (located at xq). The statistics of thi s pro cess 
is equivalent to the charge transfer statistics ( 4.16 ), as 
mentioned at the end of the previous subsection. 

The second mode, Fig. |^b, has ml = 0, mr ^= and 
wo <C mr (assuming that the scatterer is located much 
closer to the left reservoir than to the right reservoir) . It 
measures the amount of momentum transferred from the 
left to the right reservoir. Its statistics reads 



+ (e 2 £ pFMR - 
+ r(e ?PFUR 



i)[/ R -r/ R (i-/ L )] 

-1)[/l(1-/r)+/r(1 



/l)]}- (5.6) 
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I t can not be reduced to the charge transfer statistics 
(4.16) by a substitution of variables, and in particular 
does not reduce to a double Poissonian in the limit r — > 0. 
(It remains multinomial in this limit.) Comparing the 
second cumulant of momentum with the second cu- 
mulant Ch arge of charge [the terms of order £ 2 in Eqs. 
(4.16) and (|5.6[)1, we find (setting = 1) 



Time ordering then no longer reduces to input-output 
ordering and we need the Keldysh technique to make 
progress |2j| . Following the analogous formulation of the 
charge co unt ing statistics flQl , we write the generating 
function (3.4) as a single exponential of an integral along 
the Keldysh time contour, 



- { PF /efC t 



(2) 
chai 



-tp F k B T(l 

TT 



r). (5.7) 



The difference vani shes at zero temperature, in accor- 
dance with Eq. ( |5.3| ) . It is independent of the voltage (as 
long as the energy dependence of T can be ignored), so 
the diff eren ce is an equilibrium property. 

Eq. (5/7) can be given a physical interpretation by 
grouping the electrons to the right of the scattering region 
into n> right movers and n< left movers. The momen- 
tum transfer to the right reservoir is proportional to the 
sum n> + n < while the charge transfer is proportional to 
the difference n > — n < , hence 



C (2) -^l gcC x<<(->+-<) 2 }} 



4((n>n< 



(n>)(n<)). 



«(n>- n<) 2 )) 
(5.8) 



We see that the difference measures correlations between 
left and right-moving electrons. Such correlations are 
due to elect rons that are backscattered with probability 
1 — r. Eq. (5/7) describes the variance in the number of 
such backscattered electrons, given that electrons in an 
energy range k B T leave the right reservoir independently 
of each other. 



T{0 = T± exp 



, (6.1a) 



F±(r,t) = — Y] i})l{r,t)p a u a p{Y)p l 3'ip (T {r,t). (6.1b) 
m — ' 



We have written the force operator in second quantized 
form, as in Eq. ( f4.l| ) , but do not assume that the electron 
field operator ip±(r,t) = i/j(r,t±) takes its asymptotic 
form in terms of incident and outgoing states. 

The generating function can be expressed in terms of 
the Keldysh Green function G, 



|ln^) = J-J> jf dt' j dRu a0 (R) 

x ^-G aa (R,r,t',t';0 . (6.2) 

OTaOrp r=0 



VI. EVALUATION IN TERMS OF THE 
KELDYSH GREEN FUNCTION 

A scattering approach as in Sec. |^ is not possible 
if the displacement u(r) varies in the scattering region. 



The Green function G aa i is a 2 x 2 matrix in the indices 
cr,a' G {+, — } that assure the correct time ordering of 
the operators. It is defined by 



G CTCT ,(R,r,M';0 = 



-ia (r ± ^ ff (R + ±r, t)i>l (R - |r, if) exp 


|£J>' Jdr' F ± (r',t')_ 


) 


(t± exp 


\aldt'^dv' F ± {v',t') 


> 



(6.3) 



VII. APPLICATION TO A DIFFUSIVE 
CONDUCTOR 



For lon g de tection times we may assume that the Green 
function (6.3) depends only on the difference t — t — t' 
of the time arguments. A Fourier transform gives 



We apply the formalism of Sec. VI to the example of G aa i (R, p, e; £) = / dr I dr ■ 



y — ip-T — l6T 



G„ 



diffusive electron transport through a freely suspended 
disordered wire. The semi-classical calculation of the 
transverse momentum noise in this geometry was done 
in Ref. B, so we can compare results. 



■{R,r,T;0- 
(7.1) 



We write p = |p|n and use the fact that in the semi- 
classical limit the Green function is peaked as a function 
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of the absolute value |p| of the momentum. Integration 
over this variable yields the semi-classical Green function 




GW>(-R,n,e;£) = - / de' G aa ,(R,nV2me , ,e;£). 

(7.2) 

We next make the diffusion approximation, expanding 
the n-dependence in spherical harmonics, 



G GG , (R, n, e- = G^, (R, s; + n a G^ t (R, e; 



(i) 



+ (n a n p - lS a p)G^ a(Tl (R,£;0 



(7.3) 



Substituting Eq. ($L$) into Eq. (j^D we find 



d? 



InP(£) = \tE F v cr / de / dRu af3 (R) 

<T=± ' ' 



y a0 GW(R,e ] $ + ^G ( X*(R,e;O 



(7.4) 



where v = p 2 F /2ir 2 v F is the density of states. 

The equation of motion for the semi-classical Green 
function in the diffusion approximation is derived in the 
same way as for the charge statistics 10 1. We find 



2ln ^inrG+ [G (0) ,G] +& F lu a pn a n p [T 3 ,G] = 0. 

(7.5) 

The length I is the mean free path, assuming isotropic 
impurity scattering. The commutators [..,..] are taken 
with respect to the Keldysh indices cr, a' and t 3 is the 
third Pauli matrix in these indices. The Green function 
satisfies the normalization condition G 2 = I that is re- 



spected by the differential equation (7.5). The boun dary 
conditions at the left and right ends of the wire are |1(| 



G L 
Gr 



1 - 2/l 2/ l 

2-2/ l 2/ L -I 

1 - 2/r 2/r 

2-2/r 2/r-1 



(7.6) 



By projecting Eq. (7.5) onto spherical harmonics we 



find that, to leading order in l/L, the second harmonic 
G^ 2 ) depends only on the zeroth harmonic G^ : 



G ( S = \pFl{u a p - ^6 aP u yy )G^[r 3 ,G^][l + Q(l/L) 



(7.7) 



Combining this relation with Eq. (7A) we see that the 
momentum statistics of a transverse mode, with u xx — 0, 
u xy 7^ 0, follows from 

^lnF(0 = y PF lE F v J de J dR <P 

x(r 3 G^[r 3 ,G^}) . (7.8) 

V / (7(7 



It remains to compute G^-°\ To calculate In F to order 
£ 2 , that is to calculate the variance C^ 2 ' of the force noise, 
it is sufficient to know G^ for £ = 0. The solution to 
the unperturbed diffusion equation ( |7.5| ) is known ifTofl , 



GW(R, £ ;£ = 0) 



1- 2/(R,e) 2/(R,e) 

2- 2/(Re) 2/(R,e)-l 



(7.9) 



where f(R,e) = / L (e) + (x/L)[f R (e) - / L (e)]. (The co- 
ordinate x runs along the wire, from x = to i = L.) 
We find 

G (2) = t—p F lE F vA / dxdeu 2 xv (x)f{x,e)[l - /(x,e)], 
15 J 

(7.10) 

with A the cross-sectional area of the wire. This is the 
same result as in Rcf. ||. 

More complicated networks of diffusive wires, includ- 
ing tunnel barriers or point contacts, can be treated 
in the same way. In such situations the unperturbed 
Green function G(°)(R,e;£ = 0) can be determined us- 
ing Naza rov' s circuit theory p4| and then substituted 
into Eq. 



VIII. CONCLUSION 

We conclude by estimating the order of magnitude 
of the cumulants of the displacement distribution P(Q) 
of a vibrating current-carrying wire. For an oscillator 
with a large quality factor only the even order cumulants 
((Q 2k )) are appreciable, given in good approximation by 

«Q 2fc » « ^(X^o)- 2fe - lim i«AP(i) 2fc )), (8.1) 



cf. Eqs. ( |3.9D and (|3.1l| ). The cumulants of transferred 
momentum AP have been calculated for a single-channel 
conductor with a localized scatterer in Sec. M. At zero 
temperature one has 



lim h(AP(t) 2k )) = 

t— >oo t 

x 1_ ln[l + r(e £ 



d? 



2 k 



eV 
2nh 

-1)1 



P F (ur + u L 



2u ) 2k 
(8.2) 



cf. Eq. (5.3). (We have reinserted Planck's constant fi for 
clarity. ) 



Combining Eqs. (8.1) and ( p.2| ) we see that in order of 
magnitude ((Q 2k )) ^ (eV Q/huj ){p F / ' MuJo) 2k . Inserting 
parameter values (following Ref. 0)^=1 mV, Q = 10 3 , 
w /2tt = 5 GHz, p F = 5 ■ I0- 24 Ns, and M = lO" 20 kg, 
we estimate 



((Q 2k )) 1/2k « 10 4/2fc x 10" 



(8.3) 



Detectors with a 10 4 A sensitivity have been proposed 
p5[. For a measurement of higher order cumulants one 
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would want cumulants of different order to be of roughly 
the same magnitude. This can be achieved by choosing 
the number eVQ/hcuo not too large. For the parameters 
chosen above, ((Q 4 )) 1/4 /((Q 2 )) 1/2 « 0.1 . 

The theory presented in this work is more than a 
framework for the calculation of higher order cumulants 
in the momentum transfer statistics. It also provides for 
a formalism to treat quantum effects in electromechani- 
cal noise. A first application, to quantum size effects in a 
constriction, has been realized pi[ . Other applications, 
including resonant tunneling, superconductivity, and in- 
teraction effects, are envisaged. 
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APPENDIX A: DERIVATION OF THE UNITARY 
TRANSFORMATION (2.4) 



We demonstrate that t he op erator U given in Eq. (2.4) 
has the desired property (2.3) of eliminating the phonon 
displacement from the ion potential. By expanding the 
exponential in Eq. (2.4) we calculate the effect of U on 
a one-electron and onc-phonon wavefunction in the posi- 
tion space representation: 

W,g)= \\J\\ 1/2 Tp[r-qu(r),q]. (Al) 



We prove Eq. (2.3) by calculating matrix elements, 



(i/h\tfV[r - Qu(r)]U\ih) = / dv J dq || J\\ Vi> - «u(r), q]V[r - qu(r)]^ 2 [r - qu(r), q] 

= JdrJ dqil>t(i,q)V(r)ih(T,q) = (ipi\V\H- (A2) 
I 



The unitarity of U follows as the special case V = 1. 

We now justify the replacement of p = U> pU wit h p 
and A = U^AU by A in the generating function (2.8), in 
the limit of a long detection time t. Since Q commutes 
with U, it is sufficient to consider A — P. (Then A = 
A(Q, P) in the more general case that A is a function of 
both Q and P.) To first order in the displacement one 
has 



p = p - n + o(u 2 ). 



(A3) 



The difference between P and P is of the order of the to- 
tal momentum II inside the wire, which is i-independent 
in a stationary state. Since the expectation value (as well 
as higher cumulants) of P increases linearly with t, we 
can neglect the difference between P and P for large t. 

To justify the replacement of p by p we note that the 
effect of U on the initial state is to shift the electron 
coor dinates by the local phonon displacement [cf. Eq. 
(Al)]. This initial shift has only a transient effect and 
can be neglected for large t. 



APPENDIX B: EFFECTIVE MASS 
APPROXIMATION 



We start with the Hamiltonian (pTs) with V = V\ 



lat 



In the absence of any deformation of the periodic 



lattice one has, in the effective mass approximation, 

^P 2 + Mat(r) = ^P* 2 . 
2m 2m* 



(Bl) 



The quasimomentum operator p* is defined in terms of 
the Bloch function <?(r) by p* = —ig'Vg -1 . We seek 
a similar approximation to the same Hamiltonian in a 
distorted lattice, assuming that u is sufficiently smooth 
that we can neglec t de rivatives of the shear tensor u a /3- 
The Hamiltonian (2.5) (for one electron) then has the 
form 



H 



1 

2m 



Pa($af3 - 2Qu a p)p[3 + Mat + Vi v 



M 



-PU + rta"<c 



(B2) 



For small displacements Q the real symmetric matrix 
X a p = 5 a p — 2Qu a f} is positive definite. We can therefore 
factorize X = TT T , with T real. We change coordinates 
to f 



= T 1 r and find 
H 



_1 8 ^_ + Vi , Tf N 

2 m dr a dr a a 

— PU + Sla^a. 
M 



VWp(Tr) 



(B3) 



We now make the assumption of a deformation inde- 
pendent effective mass II [if], that is to say, we assume 
that the Hamiltonian with the distorted lattice potential 
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Vi a t(Tf) is approximated as in Eq. (Bl) with distorted 
Bloch functions, but the same effective mass m* . Hence 



H 



1 



2m* 



d 1 
s(Tr) — 



df a g(Tr) 



Kmp(Tf) 



- — pu + n a f a. 

1\1 



(B4) 



at the Hamiltonian 



H = ^p* a (S aP - 2Qu a0 )p} + V imp - -^Pn + fla^ 



(B5) 



Transforming back to the original coordinates we arrive given in Sec. II. 
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